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Abstract 

Using analogs of familiar image methods in electrostatics and òptics, we show how to construct 
closed form wave packet solutions of the two-dimensional free-particle Schròdinger equation in 
geometries restricted by two infinite wall barriers separated by an angle Qjy = vr/iV. As an example, 
we evaluate probability densities and expectation vàlues for a zero-momentum wave packet solution 
initially localized in a 03 = 7r/3 = 60° wedge. We review the time-development of zero-momentum 
wave packets placed near a single infinite wall barrier in an Appendix. 
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I. INTRODUCTION 



Examples of closed-form wave packet solutions of the time-dependent Schròdinger equa- 
tion (SE) in introductory quantum mechanics are at a premium. A few Standard books 
treat the case of wave packets for the harmònic oscillator (either coherent or squeezed state 
versions), often using methods which are rather advanced, such as propagator techniques 
Q]. Almost all Standard texts, however, develop the explicit example of the free-particle 
Gaussian wave packet, which requires only simple integration, and which yields a closed 
form solution. These solutions can, in turn, be used to easily evaluate many expectation 
vàlues and examine both semi-classical propagation as well as wave packet spreading. 

A simple generalization of that example to a slightly more 'interacting' case is that of the 
single infinite wall in one-dimension (or particle on the half-line) defined by the potential 



. x > 
V(x) ={ . (1) 

oo x < 



12] 



Andrews |2j has noted that simple linear combination solutions of the form 

- í ip(x,t) -ij;(-x,t) x>0 

tl){x,t) = < (2) 
10 x < 

will also satisfy the Schròdinger equation for x > 0, provided that ip(x, t) is a 1D free-particle 
solution, since the free-particle Hamiltonian is invariant under x < — > —x. This form is then 
seen to also satisfy the necessary boundary condition at the infinite wall barrier for all 
times. Such mirror solutions have been used to pedagogically examine the detailed nature 
of the 'collision' of a wave packet with such an infinite wall barrier j^J, for which many exact 
results are possible [4f, as well as motivating discussions of possible experimental tests of 
the problem of the "Deflection of quantum partides by impenetrable boundartf 

A simple generalization of this 'mirror' problem to two-dimensions involves two infinite 
wall barriers placed at right angles (90°), defined by a potential 

{0 for x > and y > 
(3) 
oo otherwise 

where the particle is constrained to be in the first (upper-right) quadrant. Product solutions, 
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using 'mirror' combinations as above, of the form 

- I [ipi(x,t) -ipi(-x,t))[ip 2 (y,t) - ip 2 (-y,t)) for x > and y > 

y{x,y;t) = < (4) 

I otherwise 

satisfy the 2D Schròdinger equation for x,y > 0, as well as the boundary conditions on 
the two intersecting infinite barriers. The form of the product wavefunction in the first 
quadrant, 

ipi(x,t)ip 2 (y,t) = i/;i(x,t)ip 2 (y,t) -^(x^^i-y^) -ipi(-x,t)ip 2 (y,t) (5) 

+ipi(-x,t)ip 2 (-y,t) 

suggests that an even more general form, namely 

4>(x, y; t) = Ç(x, y; t) - Ç(x, -y; t) - Ç(-x, y; t) + Ç(-x, -y; t) , (6) 

will satisfy the appropriate SE (for any Ç(x, y; t) which does) since the 2D free-particle 

Hamiltonian is invariant under both x « > — x and y < > —y. The additional terms in 

Eqn. (JüJ) correspond to two 'mirror' terms reflected in each infinite wall barrier (flipped in 
sign) with an additional term (with the original sign) arising from two reflections, all terms 
being necessary in order to satisfy the boundary conditions on both infinite wall barriers. 

Thought of in this way, such solutions are immediately reminiscent of image solutions 
to the related problem of two perpendicular conducting plates in electrostatics or the even 
more fundamental problem of the images of a single real object as seen in a pair of mirrors 
placed at right angles. The pattern of 'signs' used in each of those two cases is identical 
to that seen in Eqn. (JüJ). This suggests that a much wider variety of 'two- wall' problems 
in quantum mechanics can be constructed using ones intuition from such classical textbook 
examples. 

The potential in Eqn. (jïïj) is a specific case of two infinite barriers intersecting with an 
arbitrary angle 0, defined by 

ÍO ifr>OandO<#<0 
Ve(r,6)={ . (7) 

I oo otherwise 

The corresponding electrostatic (or òptics) problems for the special cases of 0jv = k/N 
can be solved by making use of a finite number of virtual 'image' charges (or true images) 
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given by the following prescription. Place the original charge (or object) at a general point 
(x,y) within the allowed wedge region. Add a 'mirror' image charge (or object) at (x, —y), 
but with opposite sign (handedness) . Then repeatedly (N — 1 times) rotate the location of 
these two charges (or objects) through 9 = 20^, keeping track of the appropriate sign (or 
handedness). For = 0^ = tt/N, this process closes on itself, and the resulting set of 2N 
charges (or objects), 2A^ — 1 of which are auxiliary ones, satisfies the boundary conditions for 
the problem. (The cases of N = 1,2 correspond to the single infinite wall and the 90° wedge 
in Eqn. (jïïj) respectively.) In the context of the 2D Schròdinger equation, such rotations of 
coordinates can be easily proved to leave the free-particle Hamiltonian invariant, so that 
a wavefunction of the form ip(x cos(0) + ysin(0), — xsin(0) + y cos(0);í) will be a solution 
if ip{x, y; t) is. This fact justifies the use of additional rotated versions of an existing free- 
particle solution in a linear combination, generalizing the result Eqn. (jüj) for the case of 
N = 2. 

As an example of this construction and its relevance to the problem of two infinite wall 
barriers, consider the case of ©3 = 60°, as shown in Fig. ^ where the prescription above 
yields the five 'image' points shown there. Using the appropriate locations, and signs, it is 
possible to construct a general solution as 

ip 60 (x,y;t) = ip(x,y;t) - 4>(x,-y;t) 

( x V3y V3x y \ / x V3y V3x y 
+ ^ {~2 " — ' T" - 2 ; *J - ^ {~2 + — ' T" + 2 ; * 

/ x V3y V3x y \ ( x V3y V3x y \ 
+ ^[-2 + —-—-ï t )-^{-2-—-— + 2> t ) ■ (8) 

This can readily be seen to satisfy the boundary conditions on the two infinite walls, namely 

ip 60 (x,y = 0;t) = = íp m (x,y = VSx;t) Wx,t . (9) 

The results for the case of G = tt/4 = 45° and others are obtained in exactly the same 
way We note that the 1D infinite square well has been discussed in terms of image methods 
Q], [3], requiring an infinite number of virtual terms, analogous to the optical case of two 
parallel mirrors. The N —>■ 00 limit of small angle wedges (Gat — * 0) is qualitatively different 
as it does not lead to bound states. 

The Standard 1D Gaussian free-particle solution with arbitrary initial central position 
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and momentum (xq and po) is given by 

lb(r\(x t) = = p ipo{x-x )/h p -iplt/2mh p -{x-x -p t/m) 2 /2t3 2 {l+it/t ) /-^q\ 

' y/^P{l+it/to) 

where 

ío = — ' Axo = 7!' and Apo = ^71 (n) 

are the spreading time and initial spreads in position and momentum (for the free 1D 
particle) respectively. 

Products of these of the form ^(G)^; ^^(G)^ ^) can then be used in Eqn. (jÏÏJ) for any 
desired initial conditions. A simple example of a spreading wavepacket, corresponding to 
(p y o,Pxo) = (0,0), initially near the còrner of a 6 3 = 60° wedge, is shown in Fig. |21 il·lus- 
trat ing the complex time-development possible. 

For the case of the single infmite wall, and the = 7r/2 = 90° wedge, such an approach 
can be continued to include the exact normalization of the wave function if Gaussian com- 
ponents are used as well as in the evaluation of a number of expectation vàlues and related 



quantities exactly 



D- 



For the cases of O = ir/N with N > 2, the normalization integrals 
involve Airy functions and so cannot be done in as compact a closed form, but experience 
with simpler cases guarantees that if the individual 'component' terms (the original and 
image contributions) are sufficiently separated in phase space, i.e. if their central vàlues of 
x ,Po are far apart in units of the fundamental position/momentum spreads, Ax and Ap , 
then the normalization factors differ from unity by exponentially small terms. 

These systems can be profitably discussed in introductory quantum mechanics courses, 
as non-trivial examples of the imposition of boundary conditions in a novel 2D geometry, 
in a free-particle context as opposed to more familiar bound state systems. In this context, 
the = 7r/2, 7r/3 and G = 7r/4 cases can be compared to bound state quantum wells with 
related geometries, as special cases of wave propagation near a 'còrner'. The N = 2 case of 
two walls at right angles provides an example of relevance to the 2D square well, where the 
very short-term evolution of wave packets in that bound-state system would be that of just 
such a 'còrner reflector', but with a more complex long-term time-dependence, including the 
possibility of quantum revivals. The 45° — 45° — 90° triangular well (or quantum billiard) 
is soluble using the odd linear combination of degenerate eigenstates of the corresponding 
2D square well, namely ip(x,y) = (u n (x)u m (y) — u m (x)u n (y)) / y2 (where m ^ n) which 
satisfies the boundary conditions on the sides of the square, but also along the y = x 



diagonal boundary. The case of an equilateral (60° — 60° — 60°) triangle infmite well is also 
known (perhaps not nearly as well as it should be) to also have simple closed form energy- 
eigenstate solutions |8] and the 3 = 60° wedge can be thought of as one 'còrner' of that 
triangular billiard. (The construction of wave packets in all three bound state systems has 
been discussed in a pedagogical context in Ref. 

The ability to manipulate relatively simple closed-form solutions enables a wide variety of 
visualizations to be easily generated, so students can examine the analogs of many classical 
'bouncing' trajectòries, including more complex ones than the simple 'in-and-out' paths for 
the N = 2 'còrner reflector'. More sophisticated projects involving the numerical evaluation 
of time-dependent expectation vàlues (to be compared with classical expectations) or of 
the momentum-space probability densities are now more within reach since a closed-form 
solution for ijj(x,y;t) already exists. For example, while a classical particle placed near the 
còrner of a wedge would remain at rest, Fig.|2]suggests that the long-term expectation vàlues 
of (x)t and (p)t will have interesting non-classical behavior. To examine just this effect, the 
time-dependent expectation vàlues (x)t and (y)t for that case are evaluated numerically and 
shown in Fig. El for the same parameter set as in Fig. 121 We note in this case that the 
long-term vàlues of (p x )t and (p y )t, evaluated using the bàsic relation (p x ) t = md(x) t /dt, 
can be seen to satisfy tan(#) = (p y )t/ (Px)t ~ 0.577 giving 9 = 30° to a very high numerical 
accuracy; thus, the wave packet develops in time in a manner consistent with its restricted 
geometry. The interplay between the generation of non-zero vàlues of the expectation vàlues 
of momenta and the conservation of kinetic energy are qüestions which arise naturally to 
students when examining such systems. These are perhaps more easily addressed in the 
context of a single infinite wall (the N = 1 discussed in Ref. 51) and we review some 

aspects of those qüestions in Appendix 1X1 

Finally, and perhaps more importantly, the notion of múltiple Gaussian wavepackets 
overlapping and demonstrating interference phenomena, such as seen in Fig. |2j is a useful 
introduction or way to motivate the study of famous experiments which have shown the 
" Observation of interference between two Bose condensates" (as first discussed in Ref. , 



but also observed in Ref. ll|, and subsequently extended to see matter wave interference 
between large numbers of BEC's Such experiments have been analyzed in terms 

of linear combinations of localized Gaussian wave packets, allowed to expand after being 
released, and exhibiting the observed interference behavior as they overlap. One can imag- 
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ine single BEC's released near one or more infmite wall boundaries (modeling atòmic 
mirrors) being described by a wave function solution much like that in Eqns. (jü)). or (jSj). 
discussed here in a more pedagogical context. 



APPENDIX A: MOMENTUM AND KINETIC ENERGY FOR A ZERO- 
MOMENTUM 1D WAVE PACKET IMPINGING ON AN INFINITE WALL 

For a classical point particle rebounding elastically from a wall, the impulsive change in 
momentum (from ±po to =FPo) and the fact that the total kinetic energy is unchanged (since 
T = T?/2m), have analogs in the interactions of a quantum wave packet with a 1D infmite 
wall (^j], 0| . The time-development of the expectation vàlues of momentum, kinetic energy, 
and the spread in momentum for a initially 'zero-momentum' wave packet expanding near 
an infmite wall can easily raise related qüestions which can be nicely visualized, understood 
intuitively, and for which long-time approximations can be derived, all aspects which we 
review in this Appendix. 

In Fig. HI we plot the position-space (\ip(x, t)\ 2 versus x) and momentum-space (\4>(j>, t)\ 2 
versus p) distributions for a simple 'mirror' solution of the form in Eqn. (J2J) using a Gaus- 
sian ip/Q){x,t) with po = 0. The momentum-space wave function is obtained by numerical 
Fourier transform of the closed-form ip(x,t). A free-particle p = Gaussian wave packet 
would simply spread, while maintaining its Gaussian form, with an increasing uncertainty in 
position given by Axt = Axo\/l + (t/t ) 2 ; the corresponding momentum-space probability 
distribution, \(f>(p, t)\ 2 , would remain unchanged, since 4>{p,t) = 4>(p,0) exp(—ip 2 t/2mh) for 
the 1D free-particle. 

For this case of the 1D infinite barrier, however, the wave packet components which move 
towards the wall (p < 0) must rebound, flipping sign and combine with the right-moving 
(p > 0) vàlues which are not flipped, yielding an overall non-zero value of the expectation 
value of momentum, as shown by the vertical dotted lines in the \<p(p, t)\ 2 plots on the right 
of Fig. 0] The negative vàlues of momentum in the initial <j)(p, 0) are, roughly speaking, 
'folded over' on top of the positive vàlues after reflection in sign, yielding a distribution 
which is therefore also narrower; a crude estimate based on this picture would be that the 
momentum-space spread would be roughly half as large as its initial value. For increasingly 
large time, any small negative components will eventually rebound from the wall so that 
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\(f>{p, t)\ 2 becomes non-vanishing only for p > 0. 

For any free-particle solution, the expectation value of the kinetic energy, given by (T) = 
(p 2 )/2m, will remain fixed, but the relation 

|p « (p 2 )o = (P 2 )t = (P) 2 t + (A Pt f (Al) 

allows for a non-zero value of (p)t to develop (as it does here, for the physical reasons 
mentioned above) but at the expense of a smaller width in momentum-space (because the 
negative components are 'flipped' on top of the corresponding positive ones.) 

The interference structure seen in the momentum-space probability density can also be 
understood from simple arguments. Momentum components labeled by p will have an as- 
sociated wavelength of À = liïhjp and as the wave packet spreads (and interferes with its 
'mirror image') there can be destructive interference between ip(x,t) and ip(—x,t) (destruc- 
tive because of the sign difference in Eqn. (J2J)) when nX = 2xq or when p = n(irh/xo). (This 
type of interference structure in the momentum-space probability density has been studied 
in the context of Bose-Einstein condensates ^| and is a general feature of the overlap of 
two (or more) expanding Gaussians.) 

An excel·lent approximation to the long-term (t » to) form of 4>(p,t) can be derived by 
first considering the linear combination of free-particle Gaussian momentum-space solutions, 
namely 



^ — a 2 p 2 /2 —ip 2 t/2mh —ipxo/h / ^ —a 2 p 2 /2 —ip 2 t/2mh +ipxo/h 



íh(r> f\ — ( > ' ' .-"-/'-/- , -'r'l ''2, 1,1, -,,,,■„,■!> i " 

2NU-^= sin (^) e -V/2 e -^t/2ma_ (A2) 



\/27r V h ) 

The two contributions differ only in their central location, determined by the exp(±ipxo/h) 
terms, corresponding to the ip(x,t) and ip(—x,t) terms. The normalization factor, N, is 
exponentially close to unity if 2x » (3 and a = (3/K. For the 'mirror' solution in Eqn. (J2J), 
the Fourier transform of íp(x,t) required to obtain 4>(p,t) initially only samples the isolated 
Gaussian term peaked at +xq, giving the Standard Gaussian form for both i/j(x, 0) and 
4>(p, 0) as shown in Fig. 0] For very long times, however, when both terms in the 'mirror 
solution' have had time to expand and overlap substantially, the Fourier integration over 
x > increasingly samples contributions from both ip(x, t) and ip(—x, t) more evenly, giving 
a result for the associated <j)(p, t) which is close to that in Eqn. (jA2j) . but with an additional 
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factor of \/2 (to ensure proper normalization) since it's eventually restricted to p > 0. This 
gives an approximate expression for the momentum-space probability density of 



\4>( P ,t »í )P 



4a fpxo\ z 

^F sm hrJ e 



(A3) 



for p > and zero otherwise, and this form exhibits zeros at integral múltiples of p = ttK/xq. 
This limiting case is also plotted in Fig. 0] (as the dashed curves on the right) and the 
numerically obtained results do approach that form for long times. The long-term value of 
the expectation value of momentum can also be approximated by noting that 



(p) 



t»t 



+oo 



V 



4a . fpx \ 2 _ a 2 2 
sm — — e p 

'7T \ k J 



dp 



(A4) 



where we approximate the oscillating úv?{pxo/h) term by its average of 1/2. This then gives 



í»ío 



7T 



2) 



2vra 2 



Tï 



7Ï 



h 



V2(3 



O.QAp 



(A5) 



which is what is observed numerically for t » to and which is close to the crude estimate 
mentioned above. (Note that a similar expression was found in Ref. [3] for the spread in 
position of a Gaussian wave packet as it hits an infmite wall and is temporarily compressed.) 
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FIG. 1: Image construction for the 03 = tt/3 = 60° wedge potential, leading to the solution 
Eqn. ©. 
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(b) 




FIG. 2: Contour plot of the probability density (\ip(x, y; t)\ 2 versus (x,y)) for a spreading wave 
packet in the O3 = 60° wedge, using the solution in Eqn. JHJ) with individual terms being a 
product of Gaussian forms as in Eqn. (|10j) . The example shown here has (p x o,Pyo) = (0,0) and 
(xo,yo) = (5,3). Numerical vàlues of h, m, /3 = 1 are used which give £q = 1 as well. The (a)-(d) 
cases shown correspond to t = 0, 5, 10, 15 in these units. 
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FIG. 3: Expectation vàlues, (x)t and (y)t versus t, for the zero-momentum initial state shown in 
Fig.[21 The insert shows the short-term time-development, initially consistent with the (stationary) 
expansion of a zero-momentum wave packet, until the higher momentum components begin to 
reflect from the infinite wall boundaries, giving non-zero vàlues of (p x )t an d {Py)t- 
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I\j/(x t)l 2 vs. x ^ l$(p t)l 2 vs. p 
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FIG. 4: Plots of the position-space (left) and momentum-space (right) probability densities versus 
time for a po x = wavepacket allowed to expand near an infinite wall. The vertical dotted lines 
show the time-dependent expectaion vàlues (x)t (left) and {p)t (right). The dashed curve on the 
right is the t >> íq solution from Eqn. ()A3|) . Vàlues of xq = 3 and h, (3, m = 1 are used, giving 
*o = 1- 
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